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ABSTRACT 


Climate models based on global radiative equilibrium are 
normally so complicated that they require extensive computer 
codes to provide adequate accuracy. However, by simply modifying 
the concept of a gray atmosphere, we can obtain a reasonably 
correct mean global temperature. This elementary model is then 
used to estimate the effects of changes in the abundances of 
minor infrared absorbers and changes in the solar constant or 
earth albedo. When applied to a Budyko-Sellers zonally averaged 
model, the quasi-gray model could give a physical basis for the 
latitude dependence of outgoing radiation and of opacity due to 
H 2 0 vapor content. The latter effect constitutes an important 
positive feedback on surface temperature. 
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What approximate, analytic models lose in accuracy, compared 
with more elaborate computer-based analyses, they gain in clarity 
and adaptability. The model developed here is a case in point. 

It is a simple matter to change the CC^ content of the atmosphere 
and f : nd what change will occur in the global mean temperature. 

One can go a step further and find how the H^O vapor content will 
be altered and include that feedback in the solution. An extension 
of this technique, of approximating the atmosphere as a quasi-gray 
one, could improve the utility of analytic, zonally averaged 
models, such as North's (1975) version of the Budyko (1969) and 
Sellers (1969) models. 

The troposphere is not in radiative equilibrium nor is it 
even approximately gray (meaning that the absorption coefficient 
is independent of frequency) . Further, it is not vertically 
homogeneous, since H^O is distributed with a scale height of 
about 2 km, CO 2 with 8 km, and 0^ existing mainly in the strato- 
sphere. Hence it is surprising that a gray, radiating homogeneous 
model has any value at all. Certainly this kind of model would 
not yield a very useful temperature-height profile. Such a 
simplified model can be useful, however, for exploring the impor- 
tance of changes in the greenhouse effect, caused by changes in 
the atmospheric composition, and of changes in the incident solar 


flux. 
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1. The Gray Model 

In the gray model only two temperatures are important: The 

effective emission temperature, T , of the planet is fixed by a 
balance between the incident solar flux and the planetary thermal 
emission. For a rotating planet of radius R 

a o 

4ttR OT = (1 - A) ttR (tt9) (1) 

e 


where A is the omnidirectional albedo and rr ^ is the solar flux 

outside the atmosphere. The surface temperature, T , which is 

s 

assumed to be the immediate source of radiation as far as the 
atmospheric heating is concerned, is related to T e in a gray 
model by (Chamberlain, 1978, p 11) 


4 4 3 

T = T (1 + 7 t ) 
s e ' 4 g 


( 2 ) 


where r is the vertical opacity of the atmosphere, which we now estimate. 
9 

If the ground emits isotropically with intensity B (T ) , the 

A. S 

flux transmission directly to cold space is 


aT Q(t) = f d\ B (T ) J d0 P e 
S 0 X 0 ' 0 


~t/h 


4 d u 


aT 2e_ (t) 

S 3 


4 - Rt 

% J T e 
s 


(3) 


where E^(t) is the exponential integral, which we approximate with 
a simple exponential. As t varies from ® to 0 and the transmission 
function Q = 2E 3 (t) varies from 0 to 1, the coefficient 4 varies 
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from 1.2 to 2.0. For H 2 <0 and C> 3 we choose 9 = 1.66, which gives 
the exact value of E 3 when Q = 1/2; for C0 2 at 15 jLtm, we take 
S * 4* s 1.2. 


For a realistic atmosphere the same kind of considerations 


hold except that r is a function of and 

n 


Q « 


Pi e 


-ftr 


(4) 


j=l 


where the summation is over narrow wavelength intervals and p^ is 

til 

the fractional flux in the i — interval. 


Pi = 


TTB li (T s )AX i 

CTT 4 


(5) 


To construct realistic temperature profiles, even for a 

truly gray atmosphere, one has to test the temperature gradient 

for adiabatic stability. In this investigation, that complication 

does not enter, since it would not change T g and would have only 

a small effect on T . 

s 


2 . Earth's Infrared Absorption 

Between 0.54 and 8 j^m has strong vibration-rotation bands 
and longward of . 1 5 occur rotational lines from the ground 
vibrational state. The Earth's thermal radiation is mainly 
confined to the 5 to 30 um region. The H^O molecule is thus a strong 
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absorber over the thermal spectrum, except for the important 
8 to 18 ym interval, where the Earth's emission peaks (at 290 °K, 
*max = ^ * Tlle to M* 11 is dominated by the CC> 2 

v 2 ~ (bending) mode fundamental, leaving an 8 to 12 um transparency 
window, a portion of this window is blocked by the 9.6 jum band of 
0 2 - Greenhouse heating of the atmosphere is thus dependent on 
composition in several distinct ways. 

First, water vapor is dominant over a wide wavelength interval, 
although in much of this region it is not optically thick. Hence 
atmospheric cooling to space is critically dependent on H 2 0 
abundance and even the vapor in the stratosphere is significant:. 

(The latitudinal variation of H 2 <0 is the main cause of the dif- 
ference in temperature profiles between polar and tropical zones.) 

Second, C0 2 fills a large part of the H 2 0 window, and current 
C0 2 abundances are enough to make the 15 ym region optically 
thick. Major CC> 2 abundance changes will alter the greenhouse 
heating mainly by CC> 2 absorption in the 10 ym bands. 

t 

Finally, minor constituents that absorb strongly in the 
8 - 12 ym region can be critically important. Nitrous oxide (N 2 0) 
and a wide array of hydrocarbons are potentially important, and 
the gray model can oe used to estimate their effects. 

The abundance of water vapor typically varies from 0.2 to 
2 

2 gm/cm (or cm of precipitable water) . We adopt a global average 
2 

of 1 gm/cm . With sufficient accuracy H 2 o may be regarded as 
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corapletely opaque except in the near infrared and visual 
(p^ ss 0.03); from 8.5 to 12.5 where r 0.10 and p = 0.29; and 
12.5 - 15 /Ltm where r < 1 and p = 0.14. However, in the latter 
region H^O does not contribute to atmospheric heating since it is 
dominated by C0 2> and the ^0 therefore can be treated as trans- 
parent from 12.5 to 18 j^m with p = 0.26 (Kondratyev, 1969, p 119). 
Thus Eq. (4) gives the absorption function due to water alone as 

A(H 2 0) = 1 - Q(H 2 0) = 0.45 (6) 

As C0 2 increases, the mean opacity within a fixed wavelength 
interval also increases, but also the band effectively widens as 
weaker absorptions become stronger. Here we will not consider the 
broadening of the 15 /jjm band because a more important (optically 
thin) effect is the emission by the 10 /urn bands. Thus the treat- 
ment of 15 fjjm itself is not quite correct, but the treatment of 
15 fim + 10 /im together should give a good representation of CC> 2 . 

For C0 2 a volume mixing ratio of f = 333 ppm corresponds to 
an equivalent thickness at STP conditions of 

§ (C0 2 ) = 266 atm-cm (7) 

An empirical formula for the mean transmission in the 
nearly saturated 15 n m band is (Kondratyev, 1969, p 126) 

^Q(15 M )') = exp (-S*a-k) 


( 8 ) 
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where the absorption data near § = 300 atm-cm yield /?* = 1.2, 
a * 0.361, and b * 0.23. Then the total absorption due to C0 2 
at 15 fim (or more precisely between 12.5 and 18.2 jim) is 

A(C0 2 ) = 1-Q(C0 2 ) = p(15 M )[l - exp(-£*a* b )] = 0.27(0.791) = 0.214 

(9) 

For the weak CC> 2 band absorption we use the treatment in 
Section 4, below, for minor constituents. The properties of the 
two intercombination bands at 10 fjjn are listed in Table 1, and 
from Eq. (24) we have A (10^) = 0.025. 

For % =0.3 atm-cm of ozone in the stratosphere we find 
(Kondratyev, 1969, p 137) over the band (9.4 to 9.9 /im) a mean 
transmission of 

^Q(9.6 n)) = exp(-2.8 = ) (10) 

or 

A(0 3 ) = 1 - Q(0 3 ) = p(9.6y) (1 - e“ 2,8? ) = 0.035(.57) = 0.020 

(ID 

The computation of the various A's is illustrated in Fig. 1. 

If the various substances do not overlap seriously in their 

spectral regions of absorption, the A's can be added. Thus from 

and (24) 

Eqs . (4) , (6) , (9) , -em* (11)/ we have (with 8 = 1.66) 

-£t 

Q net = e 9 = 1 - [A(H 2 0) + A(C0 2 ,15 jx) + A(C0 2 ,10y) + A(0 3 )] = 0.29 


( 12 ) 
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giving 

T * 0.748 (13) 

9 

With an albedo A = 29%, Eq. (1) gives T = 257 °K. Then 

e 

Eqs. (2) and (13) yield 

T = 287.2 (14) 

s 

which is close to the world average (~ 290 °K) . It is probably 
fortuitous that the agreement is so good, although the agreement 
does suggest that the gray model may provide some insight into 
climate changes due to changes in r . 


3 . Variation of Ground Temperature with CO ^ Abundance 
Differentiating Eq. (12) gives with Eqs. (9) and (24), below, 
dA (CO_ ) 

^ — - p(15 4 >exp(-(J*a; b ) /S*ab ? b (± ^ ) 

+ 1.08 x 10‘ 3 S ^ (15) 

V 0 dt 


In terms of the mixing ratio or volume concentration 


f £r _ rco I _ pilSj t ) j? *ab ex p(-y3*ag b ) ? b 

df ' lCO 2 ] d[C0 2 l ' 


» 4.21 x 10" 2 + 7.5 x 10 -2 = 


1.08xl0“ 3 S N* 
y 0 



0.117 


(16) 


Similarly from Eq. (2) we obtain 



9 


dT 


s 


dT 

and therefore 


3T 


16(1 + .75r ) 
9 


3/4 


3T 

s 

16(1 + . 75t ) 

g 


= 34.5 


(17) 


dT 

f j-jT - 4.0 °K (18) 

.if the CO 2 fractional concentration were changed by a small amount. 
If the CC >2 content were doubled, the gray model would give 

T = 290.0 or AT = 2.8 °K. This value compares nicely to values 

obtained from more elaborate models. In reviewing a number of 

recent studies, Schneider (1975) concluded that a doubling of 

CO. ould produce a AT of 1.5 to 3 °K. As CO- is increased 

furtner, the temperature will continue to rise; although the 

15 jum bands are already nearly saturated, the 10 urn bands are 

effectively closing the 8 - 12 ^m window, as illustrated in Fig. 2 . 


4. Radiative Effects of Minor Constituents 

The addition to the atmosphere of a minor constituent that 

absorbs in the 8 to 12 ^m window could be important. In this region 
4 

ttB /aT as .060 - .070 //i and hence Eq. (5) is 

A. S 


p. as 0.065 AX 
* M 


( 19 ) 


where AX is the width of the absorption in microns. We will 
estimate the minimum amount of a trace gas that could produce a 
1 °K change in global temperature. 

For an optically thin amount of absorber, we have simply 
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Table 1 

ABSORBERS IN 8 - 12 WINDOW* 


Species 

Band, X(am) 

f (ppm) 

Nitrous oxide 

7.8 

0.28 

CH. Methane 
4 

v . 7.66 
4 

1.6 

NH^ Ammonia 

v 2 10.5 

6 x 10 -3 

CF 2 C1 2 F-12 

v 1 9.1 

1 x 10 -4 


v 6 8.7 



v 8 10.9 


cfci 3 F-ll 

v L 9.2 

1 x 10 -4 


v 4 11.8 


Hydrocarbons 

— 

<4 x 10~ 3 

C0 2 

2v 2 ~v 3 9.4 
v^-v^ 10.4 

333 


10 


S (cm) 
-17 

-18 


7 x 10 
~ 10 _2 °(?) 


5 x 10 


-17 


3 x 10 


-17 


5 x 10 


-17 


3 x 10 


-17 


6 x 10 


-17 


<5 x 10 


-17 


8.6 x 10 


2.7 x 10 


-22 


-21 


* For uniform mixing the integrated abundance is oA/ = f cV 

atm 

25 -2 

= f x 2.15 x 10 cm Since the STP thickness is % = /N^ 

= f /2.687 x 10 19 , we have 

atm 


f{ppm) * 1.25 % (atm-cm) . 
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”^ T i 

e j * #>i<V 


( 20 ) 


where <t> is the mean value in the band, or 


<T> " U J c ax 

which gives 

A = ( .065) <*) J a dX 

= .108 S 

K 


( 21 ) 


(22) 


2 

where is the integrated column density (molecule/cm ) and S 

K 

the integrated absorption coefficient or strength (in units of 
2 

cm micron) . The more common units for strengths are cm (cross 

2 -1 
section in cm ; spectrum in cm ) : 


S = I* a d> 

\J J 


(23) 


Hence at 10 ^m, 


A = 0.108 ^ S X 2 10" 4 

v M 

= 1.08 x 10” 3 gtf S 


(24) 


_i ? 

The maximum band strengths are not likely to exceed 5 x 10 cm 


From Eq. (12) we have, for A = 1 - Q 

net 


a . - 

dT g 


= 1.66 (0.25) = 0.415 


( 25 ) 


and with Eq. (17) the T - A relation becomes 

s 
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dT dT /dr 

ST * A/dr l * 83 • <26> 

Thus for 1 K changes in T , we are searching for increments in A 

S 

-2 

the order of 1 x 10 , or from Eq. (24), 

A*/ > L-XJ£ — = 2 x 10 17 cm' 2 (27) 

~ 1.08 x 10 x 5 x 10 

25 -2 

The integrated column density of the atmosphere is 2.15 x 10 cm 
so we are concerned with mixing ratios of minor constituents 
greater than 


Af ~ 10 x 10 


(28) 


or a few parts per billion (ppb) for substances that have a single 

strong absorbing band in the 8 to 12 ^m window. 

The halogenated methanes (Freons) have been cited as a 

striking example of the effect we are examining here (Ramanathan, 

1975) . The two Freons F-ll and F-12 together have five bands 

2. *7 

with strengths ~ 3 x 10 cm. The 1975 abundance was 0.1 ppb; 

2 

an increase of 10 times would place these substances in the climate 
critical category. 

Most substances will not have bands in the window with 
strengths close to the maximum. From Eq. (27) the pertinent 
quantity for bands in the window is 

A<A/* s > 10 or Af(ppm)*S > 5 x 10 ^ (29) 

r<*> ^ 
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if an absorber is to produce a 1 °C change on future climate 
(see Table 1) . Both N 2 0 and CH 4 sit in regions where there is 
already strong absorption and they have to be analyzed (as we 
did C0 2 ) allowing for the present near-saturation. 

All of the substances in the Table except the hydrocarbons 
have been examined with a more elaborate radiative model (Wang, et al., 
1976) and found to be marginally important, if their abundances 
increase. For N 2 0, CH^, and NH^# factors of two are important. 

More attention needs to be directed toward the global 
abundance and long term growth rate of hydrocarbons, which are 
among the primary urban pollutants, and whose emission is largely 
uncontrolled. The abundance in Table 1 applies to the total 
hydrocarbon population; the strength listed is the maximum 
allowable and probably is much smaller. Aldehydes (a hydrocarbon 
family) are partly destroyed near their urban source, where they 
participate in the production of oxidants. Away from high 
concentrations of NO^., the aldehydes disappear by photodissociation 
and attachment to aerosols. 

5 . Destruction of Trace I.R, Absorbers 

Many of the hydrocarbons are not only being emitted at 
increasing rates, but their natural sinks are becoming saturated. 

The principal sink for meth- ne is hydroxyl; 
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ch 4 + OH - ch 3 + h 2 o 


(30) 


Indeed, OH plays a crucial role in scavenging a number of the less 

soluble trace gases, such as CO, H 2 S, S0 2 , and the partially 

halogenated methanes, CH Cl F and CH 3r . a depletion of OH 

x y z y. y 

may thus aggravate a build up of infrared absorbing pollutants and 
their impact on the climate. At present there is a danger of 
such a depletion occurring because of the increasing production of 
CO, which is removed mainly by 


OH + CO - C0 2 + H (31) 

9 

The present CO production from fossil fuel burning is 0.4 x 10 ton/yr, 

9 

whereas the natural production is 1 to 3 x 10 ton/yr. The 
problem is amplified by the fact that the main natural source of 
CO is the oxidation of CH^. Following reaction (30) , oxidation 
proceeds by 


CH 3 + 0 2 + m 

- CH 3 0 2 + M 

(32) 

CH 3 0 2 + NO - 

CH 3 0 + N0 2 

(33) 

ch 3 o + o 2 - 

ho 2 + ch 2 o . 

(34) 


Finally, formaldehyde is dissociated into H 2 + CO, or it may 
condense on\;o aerosols and ultimately rain out. 

The principal source of tropospheric OH is photodissociation 
of tropospheric ozone: 
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0 3 + hv(X < 4110 A) - O^D) + 0 2 (35) 

followed by 

H 2 0 + 0( 1 D) - 2 OH (36) 

Thus present levels of man-made CO emission are large enough 
to perturb the C0-0H-CH 4 cycle. Reactions that recycle OH, such 
as H0 2 + NO - OH + n0 2 ' account for some 30 percent of the 

present OH. 

6 . Feedbacks; Comments on Water and the Stratosphere 

Earth is the water planet and ’/ater affects climate differently 
in its different manifestations; vapor, clouds, sea, sea-ice, 
land-ice and snow, and foliage. Common to all models is the 
acknowledgement of the effect and inability to do anything useful 
about it. It is not unlikely that the response of water to a 
change in C0 2 is the dominant factor in climate change. 

Also physical processes in the stratosphere can feed back 
on surface temperatures. The water abundance in the stratosphere 
is only several parts per million and is probably fixed by the 
vapor pressure at the tropical tropopause [although H 2 0 can be 
made in the stratosphere by reaction (30)]. A change in the content 
of water due to a change in stratospheric temperature will affect 
the earth's greenhouse (Chamberlain, 1977). An alteration of the 
ozone abundance will change the absorption at 9.6 jim, which is 
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the main absorber in the 8 - 12 jim window. Thus an increase in 

tropospheric CH could affect stratospheric HO chemistry, the 
4 

ozone abundance, and thence the climate. 


7 . Changes in the Sclar Flux or Albedo 

Differentiating Eq. (1) and relating T to T by Eq. (2) for 

S G 

a constant opacity, we have 


AT 


1 = M(1 ~ A)tt31 

4(1 - A)tt^ 


(37) 


Consequently a fractional change of 6 percent in the solar constant 
or in 1 - A will produce a change in mean surface temperature of 


AT. 


1 5 

4 s 100 


= 0.72 5 


K 


(30) 


In practice there will be an accompanying change in H^O vapor and 
sea ice, affecting the opacity and feeding back on the albedo. 
Analysis of these effects is best handled with a latitude dependent 
(Budyko-Sellers) model. Accordingly, further discussion is 
deferred to a later paper. 
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FIGURE CAPTIONS 


Fig. 1 Thermal emission of the Earth's surface at T = 290 °K 
(after Kondratyev, 1969, p. 119), with relative values of the 
absorption, 1 - exp(-^r), indicated by shaded areas. The figure 
illustrates why trace amounts of absorbers in the 8.5 - 12.5 
region can be critically important. 

Fig. 2 Temperature versus CC >2 abundance for global radiative 
equilibrium with the quasi-gray model. 
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